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Using 230 x 10® BB events recorded with the BAR4R detector at the e+e~ storage rings PEP-II, we 
reconstruct approximately 4100 B° — > J/iI)K~^'k~ and 9930 S"*" — > J/ipK~^ decays with J/ip — »• fi^ fjT 
and e'^e~ . Prom the measured B-momentum distributions in the e'^e~ rest frame, we determine 
the mass difference m(B°) - m{B+) = (+0.33 ± 0.05 ± 0.03) MeV/c^ 



PACS numbers: 13.25.Hw, 13.40.Dk, 14.40.Nd 
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Mass (lifFcrciiccs Amj\/ = rn{M'^) — rri{M^) probe 
the size of Coulomb contributions to the quark struc- 
ture of pseudoscalax mesons M. The values of Amjvf 
for TT. K. and D mesons arc experimentally well known; 
in units of MeV/c^ they are Am^ = -4.5936 ± 0.0005, 
Atok = +3.97 ± 0.03, and Amo = -4.78 ± 0.10 [1]. 
For B mesons, Ams = (+0.37 ± 0.24) MeV/c^ [2] is less 
precise and compatible with zero. Quark-model calcula- 
tions [3] give ArriB near -1-0.3 MeV/c^ but are quite un- 
certain since the contributions from the quark-mass dif- 
ference m{d) — m{u) and from the Coulomb effects have 
similar magnitudes and opposite signs. In the case of 
Am/), the two contributions enter with the same signs. 

The value of Ams is an important input for estimat- 
ingjhe decay ratio R = r[r(4S') ^ B+B-]/T[r{4:S) 
B'^B'^] which in turn is essential for determining B^ and 
B° decay fractions at e+e" colliders where B mesons are 
produced in decays of the T(45). The leading contri- 
bution to R is given by the vector nature of the matrix 
element and by kinematics; at fixed energy it is 

Ro = [p*iB+)/p*iB')f « l + 3mBAmB/p*i , (1) 

where p*{B~^) and p*{B^) are the i?"*" and B^ momenta 
in the center-of-mass system (cms) at this energy, and 
p*Q and ms are the mean values of the two momenta 
and masses, respectively. For lAm^l below 0.5MeV/c^, 
the quark structures of T(4S') and B mesons and the 
Coulomb interaction [4] may lead to \R — Ro\ > Rq — 1. 

For measuring Ams, we use 210 fb""'^ of e+e" anni- 
hilation data recorded on the T(4S') resonance with the 
BaBar detector [5] at the SLAC e+e~ storage rings PEP- 
II [6]. Charged-particle momenta are measured by the 
tracking system consisting of a five-layer double-sided sil- 
icon vertex tracker and a 40-layer drift chamber, both 
located in a 1.5 T magnetic field of a superconducting 
solenoid. Transverse momenta px are determined with 
a resolution of about <t{pt)/pt = 0.0013 x pxc/GdV + 
0.0045 and track angles with resolutions around 0.4mrad. 

The B mesons are reconstructed in two decay modes 
with low background level: i?° ^ J/ipK'^ Tr~ and B^ — > 
J/'d)K^ [7], where J/ij) ^I'^fjr or e+e^ in both modes. 
Measurements of and J/iJj invariant masses show 
that relative momentum uncertainties Sp/p, originating 
from the limited knowledge of the magnetic field and 
the charged-particle energy losses, are below 4 x 10^^. 
A momentum uncertainty of this size leads to B-meson 
mass uncertainties of the order of 1 MeV/ . The mass 
difference Ams can be determined with much higher 
precision using _B-meson momenta because the decay 
Y{4S) BB produces B mesons with low momenta, 
p*{B) Ri 320MeV/c. At fixed cms energy y/s we have 

m^{B+y + p*^{B+) = m^iB'^y + p*^{B°) , (2) 



where Ap* = p*{B'^) —p*(B~^). The track-momentum 



uncertainties lead to 6p* < 4 x 10 



320MeV/c, 



S{Ap*) <V2 X Sp*, and, using Eq. 3, 5{AmB) < 
0.01 MeV/c^ which is two orders of magnitude smaller 
than the IMeV/c^ estimate using invariant masses. 

The energy spread of the PEP-II beams gives a i/i 
distribution with a rms width of about 5 MeV, result- 
ing in broad distributions of the true momenta Ptrue(-^) 
with rms widths of about 40 MeV/ c. The reconstructed 
p* spectra are only slightly wider since the detector res- 
olution is only a{p* — Pt^ue) ~ 15MeV/c in the selected 
B-decay modes, where a is the rms width. As input for 
Eq. 3, we use the mean values p*{B^) and p*{B^) of the 
reconstructed p* spectra. The presence of background 
prevents the two p* values from being obtained as al- 
gebraic means of all measured p* values. Instead, they 
are determined from fits with analytic functions for the 
signal and the background shapes. 

The size of a possible bias from the mean-p* method is 
estimated by Monte Carlo (MC) simulations in two steps. 
The influence of the beam smearing and the T{4S) line 
shape is studied by determining the means plj.^^(B) using 
a MC simulation in the cms without detector. We use a 
Gaussian cms-energy distribution with two parameters: 
V^mean ^^'^ ^V^' ^^'^ T{4:S) BB line shapes with four 
parameters: m(T(4S')), m{B°), m{B^), and Fq, where 
the latter is the total width at s = m'^{T{AS)). The line 
shape is parametrized following ref. [8] ; it includes initial 
state radiation, a relativistic Breit-Wigner function with 
energy-dependent width, m{B)- and s-dependent phase 
space factors, and mesoii-struc;ture effects. Because of 
the m(B) dependence of the phase space factor, the line 
shapes differ for B^ and B+. We fix Ams to either 
-1-0.3 or -1-0.4 MeV/c^ and vary the other parameters in 
the range of the results of ref. [8] . We determine p^rue (-^°) 
and p1^^^{B^) for each set of parameters and find that 
the derived Attib results from Eq. 3 are equal to the MC 
input within ±2%. The rms widths Up. of the two Pt^ue 
distributions are found to be different in agreement with 



ap.(B+)/ap.(B°)=r(S°)/r(i?+) 



(4) 



Ams = -Ap* 



p*(B") +p*(B+) 
[m(BO) + m(B+)] 



(3) 



as simple consequence of Eq. 2. 

The detector influence on the Atob bias is studied 
by a full MC simulation of generic BB decays with 
GEANT4 [9]. The simulation includes all detector and 
reconstruction effects and the same Gaussian cms-energy 
distribution as above, but uses a simpler T(4S) — > BB 
line shape with fixed F(s), without initial state radiation, 
and without meson-structure effects. The results on the 
means of p* — p*rue given in the discussion of the 
systematic uncertainties. 

The same GEANT4-based MC simulation is used to 
determine the selection criteria for B reconstruction and 
to find the fit-function types for signal and background 
in the reconstructed p* spectra. The J/ip decays into jiji 
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FIG. 1: The p* distributions of the selected data (dots with error bars) and of the selected MC events (histograms stacked on 
top of each other). The left side is for B+, the right side for the top for J/ij) — » /xyu, the bottom for J/i/; — > ee. 



and ee arc studied separately in order to control the influ- 
ences of bremsstrahlung in the ee channel, simulated by 
PHOTOS [10] and GEANT4. Muons are identified using 
a neural network with a high efliiciency of 0.90 accepting a 
rather high probability for pion misidentification (misid) 
of 0.08, while electrons are identified using a likelihood 
selector with an efficiency of 0.95 and a pion-misid prob- 
ability of 10""^. Electron tracks are combined with up 
to three nearby photons into electron candidates using a 
bremsstrahlung-recovery algorithm. Pairs of electrons or 
muons with opposite charge are fitted to a common ver- 
tex. All pairs with a vertex fit probability P > 10~* and 
an invariant mass between 3.057 and 3.137GeV/c^ are 
selected as J/^p candidates. Because of background from 
two pions in jet-like e+e~ — > qq events, we also require 
\cos{6h)\ < 0.9 for J/t/j jiji candidates, where Oh is 
the angle between one muon and the B candidate in the 
J/i/' rest frame. Since the pion misid is much lower for 
electrons, this cut is not applied in the J/i) ee mode. 
In the 5° mode with J/^ ^fjt we require in addition 
that the normalized second Fox- Wolfram moment i?2 [11] 
of the event is less than 0.4. 

Charged kaons are identified using a hkelihood selec- 
tor, based on the DIRC system [12] of BaBar, with an 
efficiency of 0.95 and a pion-misid probability of 0.05. 
The Kit pairs are formed from two oppositely charged 
tracks, one identified as a kaon and the other as a pion; 
the fit to a common vertex must give a fit probability 
P > 10~^. For suppressing background, we require an 



invariant mass m(fr7r) = m(jr*°) ± 75MeV/c2. The 
and B^ candidates are formed by combining the J/^jJ with 
the Kir-pskU candidates and with charged tracks identi- 
fied as kaons, respectively. We also require a fit probabil- 
ity P > 10~* for the common vertex. The B candidates 
are further selected by their value of AE* = E*^ — \/s/'2,, 
where E*^ is the energy of the B candidate in the cms. 

To optimize signal versus background in the p* distri- 
butions and to account for bremsstrahlung, we have cho- 
sen four different Ai?* selection criteria. For J/ij) jifj, 
we choose |Ai;*| < 55MeV for the B+ and |Ai;*| < 
25 MeV for the For J/il) ee we take -60 < l^E* < 
50 MeV and -30 < A£;* < 20 MeV in B+ and S° decays 
respectively. For the B^ this corresponds to ±3 rms 
widths of the signal, for the B^ to ±1.5 rms. The tighter 
criteria in B^ decays, where the background is an impor- 
tant contribution to the final systematic uncertainty on 
Ams, axe justified by the negligible correlations between 
AE* and p* and by the MC validation as described be- 
low. After applying the AE* criteria to the B candidates, 
there remain events with more than one candidate. The 
fraction is negligible for i?"*" (0.10% of all events) but is 
1.5% for B^. If there are multiple B candidates in the 
event, we choose the one with the best B-vertex fit. The 
selection criteria for data and MC events are identical 
with one exception: In the data, because of a bias in 
the J/ij) mass reconstruction owing to track-momentum 
uncertainties, the lower and upper limits for m(/i/i) and 
m(ee) are shifted by — 2MeV/c^. 
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Figure 1 shows the p* distributions of the selected 
data and those from the MC simulation. The MC dis- 
tributions are normalized to the data between 0.12 and 
0.45 GeV/c. They contain contributions from four classes, 
Class 1, "pure signal", candidates where all tracks origi- 
nate from true B-decay particles into the given mode 
and where the decays contain no photons includ- 
ing those combined into electron candidates with the 
bremsstrahlung-recovery algorithm, 
Class 2, "signal with radiation" , like pure signal, but 
with at least one photon from bremsstrahlung gen- 
erated bj^PHOTOS or GEANT4, 
Class 3, "BB background", candidates from B decays 

other than from classes 1 or 2, and 
Class 4, "g^ background", candidates from non-SB 
events. 

The third class also contains some signal events with 
wrong matching of reconstructed and generated tracks. 
As can clearly be seen, the BB background in 5° de- 
cays is larger than in iJ+ decays and the fraction of can- 
didates with bremsstrahlung is larger in the J/tp ee 
than in the J/ip /i/i mode. Note that, in spite of ob- 
served differences in the invariant lepton-pair mass and 
in the AE* distributions for ee and /x/i, there is almost 
no difference in the shape of the p* distributions. In Fig- 
ure 1, differences between data and simulation are seen 
on both edges of the signal peaks. They may arise from 
imperfections in describing the beam energy spread and 
the T(4S') line shape which influence 5° and iJ+ decays 
equally; the following data analysis has to account for 
the imperfections. 

The mean values of the four p* (B) spectra are obtained 
from fits. The form of the fit functions is obtained from 
the MC spectra for "pure signal" and the sum of BB and 
qq backgrounds separately. For the signal, we flnd that 
a double-Gaussian function S{p*) with six parameters is 
adequate. Its parameters arc: the number N of signal 
events (sum of classes 1 and 2), the mean p* and the rms 
width o-p. of S{p*), the fraction / of the subdominant 
Gaussian function, the peak-position difference A and 
the width ratio To- of the two Gaussian functions. The 
fits of S{p*) to the "pure signal" contributions in the 
four spectra of Figure 1 are of good quality. The fit- 
parameter values axe similar in all four spectra; only the 
(Tp. values are slightly larger in ee than in /x/x decays. It 
has been checked that S{p*) with the same parameters 
as for "pure signal" also describes the p* distributions of 
"signal with radiation" for both B^ and decays. 

The backgrounds for B'^ and B~^ are very different, re- 
quiring two different function types Uo{p*) and U-^{p*). 
We find that polynomials are adequate in both cases, lin- 
ear for 5+ and of fifth degree for 5°. The polynomials 
are determined by fits to the sum of the MC background 
histograms. Because of the complication with the mis- 
matched signal MC events, we have to use fit functions 
Uo,+ {p*) + S{p*) for determining the background poly- 



TABLE I: Results for fitting the sum of signal and background 
functions to the four data p* spectra in Figure 2. The param- 
eter N is the number of signal events and r^g is the ratio of 
the observed to the simulated background level. The values 
for p* , Op" , and A are in MeV/ c. 

B°, nij. B+, B°, ee B+, ee 

N 2280 ±50 5580 ±70 1820 ±40 4350 ±70 

p* 316.8±0.9 321.6±0.6 314.7±1.1 321.1±0.7 

o-p. 43.0 ±0.8 44.4 ±0.5 44.3 ±0.9 45.4 ±0.6 

/ 0.79 ±0.04 0.78 ±0.06 

A -51±7 -48±10 

ro- 1.46 ±0.08 1.48 ±0.08 

rbg 1.16 ±0.10 1.01 ±0.11 1.08 ±0.11 1.93 ±0.24 



nomials, where S{p*) is the best-fit signal function with 
free normalization. 

In the fits of S{p*) + C/o,+ (p*) to the p* distribu- 
tions of real data, we choose binned maximum-likelihood 
fits between 0.12 and 0.95GeV/c with bin widths of 
0.015 GcV/c. The background polynomials arc used with 
free normalizations rbg but with shape parameters as 
given by the MC fits. In the signal function, all six pa- 
rameters are left free since the signal shapes differ in data 
and MC because of the imperfect MC simulation. Since 
the \/i spectrum dominates the shapes of the two p* 
spectra, the parameters /, A, and are constrained to 
be equal for B° and B+. 

Before fitting the real data, we apply the fit to p* dis- 
tributions of the MC simulation. We divide the sample 
of reconstructed MC events in five parts of equal size, 
each with the same integrated luminosity as the data. 
The 10 fit results, combining J/ip — > /i^ and ee, have a 
mean of Ap* = p*{B°) - p*{B+) = -4.7MeV/c with a 
rms of 0.4MeV/c in good agreement with the MC input 
of -5.1McV/c. 

Figure 2 shows the p* distributions of the selected 
data events together with the best-fit functions. The 
fit results are given in Table I, and the derived Ap* 
values are (—4.8 ± l.l)MeV/c for the J/-^ jiji and 
(-6.4 ± 1.3) MeV/c for the J/-^ ^ ee mode. The two val- 
ues are consistent, we therefore use the weighted mean 

Ap* = (-5.5 ± 0.8) MeV/c 

for the final result. Before converting Ap* into Ams, we 

present a number of cross-checks and the estimates of all 
contributions to the systematic uncertainty. 

The Ap* results from different run periods of the ex- 
periment arc in agreement with each other and no charge 
dependence is observed. We find p*{B^) — p*{B~) = 
-0.3 ± 1.3 (-1.0 ± 1.4) MeV/c for the n/j, (ee) mode and 
p*{B^) - p*{B^) =0.4 ±1.8 (-1.2 ± 2.1)MeV/c for ^/x 
(ee). Varying the AE* requirements for the B candidates 
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FIG. 2: The fitted p* distributions in data, left side for right side for top for J/ip nn, bottom for J/-^ — » ee. 
The dashed hnes show the background polynomials; the dotted lines for S° show the changed backgrounds for the systematic- 
uncertainty estimate. 



by factors of 1.4 up or down changes the central value of 

the Ap* result by less than half a standard deviation. No 
sizable effect on the central value is seen when removing 
the requirement on the muon angle Oh for the J/tp can- 
didates, on m{K7r), or on the event-shape parameter R2. 

The contributions to the systematic uncertainty of Ap* 
axe summarized in Table II. The influence of the chosen 
parametrization for the signal fit-function is estimated 
by using modified parametrizations. First, we allow /, 
A, and To- to be different for and which re- 
sults in Ap* = (-5.4 ± 0.8) MeV/c. Second, we use 
one parameter less than in the nominal fit requiring 
ap-.{B°) = (Tp.(B+) X p*{B+)/p*{B^) from Eq. 4 result- 
ing in (—5.7 ± 0.8) MeV/c. We use the observed average 
variation of the three fit- method results in data and in the 
five MC validation subsamples as an estimate of the sys- 
tematic uncertainty for the signal fit-function. Since the 
backgrounds are small, we also determine p* as algebraic 
means of the four p* spectra between 0.12 and 0.45 GeV/ c 
after subtracting the best-fit background functions. The 
results agree with those in Table I within ±0.1MeV/c 
except for p*{B^\ ee) where it is 0.3MeV/c lower. 

The influence of the background in the p*{B^) spec- 
trum requires special care and was investigated by three 
methods. First, we compare the fit results for various 
AE* cuts with rbg free and rbg = 1- Second, in order 
to control the influence of a slightly different background 
shape in the signal region, we flt the B'^ data using mod- 



ifled functions Uo{p*) with the arbitrary shapes of the 
two dotted lines in Figure 2. Third, we select wrong- 
sign candidates in the channel J/'ijjK^-K^ with all selec- 
tion criteria as for the nominal candidates includ- 
ing those for 'm{KTT). The ratio Q of sclccitcci data and 
MC events is well approximated by the linear function 
Q = 0.30 + 0.78 X p*c/GeV. The function Q x Ua{p*) 
is then fltted to the selected B^ data with rbg floated. 
The second and third method give comparable shifts in 
p*{B^) and we take them as systematic uncertainty for 
the background- function; the shift in the flrst method is 3 
times smaller. Variations of the fit binning from the nom- 
inal 15 MeV/c width to 5, 10, and 20 MeV/c have a negligi- 
ble influence. The transformation from laboratory-frame 
momenta to cms momenta has negligible influence, even 
by varying the applied boost by the flve-fold of its rms 
in PEP-II. The detector influence on the Am^ bias is 
estimated by using the MC results for the means 6 of 
p* — Pti-uc as estimators for the uncertainties. The results 
are S{B")-5{B+) = (-0.14±0.13) MeV/c for the ;U;U and 
(0.25 ± 0.19) MeV/c for the ee mode. We conservatively 
use the sums of central value and rms of these results for 
the last-line entry in Table II. 

Adding all systematic uncertainties in quadrature and 
taking the larger of the two estimates (ee) leads to 

Ap* = (-5.5 ± 0.8 ± 0.5) MeV/c . 

Inserted into Eq. 3 and using 319 MeV/c and 5279 MeV/c^ 



TABLE II: Summary of the systematic uncertainties for the 
measurement of Ap* in MeV/c. 





MM 


ee 


Signal Fit- Function 


0.12 


0.17 


Background Fit-Function 


0.01 


0.03 


B° Background Fit-Function 


0.25 


0.16 


Histogram Binning 


0.08 


0.08 


Detector Bias 


0.27 


0.44 


Quadratic Sum 


0.40 


0.51 



for the mean values of B momentum and mass, we obtain 

AruB = {+0.33 ± 0.05 ± 0.03) MeV/c^ . (5) 

Contributions to the systematic uncertainty (in MeV/c^) 
come from Ap* (±0.031), the track-momentum uncer- 
tainty (±0.011), and the mean-p*-method bias (±0.007). 
The contributions from the uncertainties on the T(4S') 
boost and the B-meson mass are negligible. 

The AniB result in Eq. 5 is compatible with the present 
world average [2] of (0.37 ± 0.24) MeV/c^ but the error is 
a factor of 4 smaller. The significance of Ams being 
positive exceeds the 5a level. Inserting our Ams result 
into Eq. 1, we obtain Rq = 1.051 ± 0.009. The measured 
value of R is 1.037 ± 0.028 [2]. Given the agreement 
between these two results, we do not observe significant 
Coulomb or quark-structure contributions [4] to R. 
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